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Abstract: We use more than two years of data from the Pierre Auger Observatory to search for 
anisotropies on large scales in different energy windows. We account for various systematics in the 
acceptance, in particular due to the array growth and weather variations. We present the results of anal- 
yses and consistency checks looking for patterns in the right ascension modulation of the cosmic ray 
distribution. No significant anisotropies of this kind are observed. 
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Introduction 

The Pierre Auger Observatory has accumulated an 
important statistics of cosmic rays (CRs) above 
10^* eV = 1 EeV. The use of water-Cherenkov 
tanks allows to observe a large fraction of the sky. 
Upper limits were set on the possible fluxes from 
localized regions near the galactic center [1]. It 
is not known yet whether the sources of cosmic 
rays with an energy just below the ankle are galac- 
tic, or whether they have a cosmological distribu- 
tion [2, 3]. If the origin of cosmic rays below the 
ankle region is galactic, their diffusive escape from 
the galaxy may be efficient enough so fliat the sky 
distribution of 10^^ eV cosmic rays is not com- 
pletely isotropised as seen from the Earth. The 
predictions for the shape and amplitude of the cor- 
responding anisotropy are very model-dependent, 
but a %-level modulation is conceivable. If the 
origin of cosmic rays is extragalactic, then their 
sources should be cosmologically distributed for 
all energies below the GZK attenuation, and no 
large-scale pattern should be detectable except per- 
haps for a CMB-like dipole with an amplitude of 
0.6%. 

The AGASA collaboration found a ^ 4% ampli- 
tude modulation in right ascension (RA) in the spe- 
cific energy range 1 < E < 2 EeV [4]. The 
statistics accumulated by the Auger Observatory 
allows to study %-level large-scale patterns, but 
this is challenging due to the difficulty of control- 
ling the sky exposure of the detector with the cor- 



responding accuracy. Indeed various acceptance 
effects, such as detector instabilities and weather 
modulations, can generate spurious modulations 
in the exposure. We present here first results of 
complementary analyses in the EeV energy range, 
which set upper limits to a possible large-scale pat- 
tern. We also search for large anisotropies above 
10 EeV, where it is widely believed that the CR 
sources are extragalactic. All the analyses pre- 
sented here use the high-level trigger (T5) events 
with a zenith angle 6 < 60°, and the periods of 
unstable acquisition are removed. 

Fourier time analysis 

A modulation of the event rate over right ascen- 
sion or sidereal time may be due to a true sky pat- 
tern, but also to the combined effect of diurnal and 
seasonal variations of the detector acceptance. To 
disentangle these effects, a widely used method is 
to compare sidereal modulations to solar and anti- 
sidereal ones [5]. In the following analysis, we use 
the events recorded at the Auger Observatory from 
January 2004 to March 2007. For each frequency 
we compute the associated Fourier amplitude using 
the distribution of the times ti of the events. Rather 
than using directly ti we apply the analysis to the 
modified times imod(*) = ti + RA(z) — RAo(i), 
where RA{i) is the right ascension of an event and 
RAo{i) is the associated local sidereal time [6]. 
Using such a definition, the sidereal modulation of 
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Figure 1: Fourier analysis of the modified times 
tmod for all events from 2004 to March 2007 
(-E-mcdian ~ 0.6 EeV). Frequencies are in unit of 
the sidereal one. The dotted lines represent, from 
left to right, the anti-sidereal, solar, and sidereal 
periods. 
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Figure 2: Amplitudes (circles) and corresponding 
95% upper limits (triangles) of the first harmonic 
analysis in sidereal time vs. the energy threshold 
using the East- West method. 



^mod is equal to the RA modulation of the event 
rate. 

The Fourier analysis on all the data (5 x 10^ events 
with a median energy 6 x 10^^ eV), presented 
in Fig. 1, does not reveal a specific modulation 
at the sidereal frequency. A 3.2% solar modula- 
tion of tmod is observed, and can be interpreted 
as due to weather effects [7] and to the array de- 
ployment and maintenance which take place dur- 
ing the day. Defining Tcxp as the duration of the 
experiment, the frequency resolution is ^ l/T^xp, 
which is enough to resolve the solar and sidereal 
frequencies (as can be seen from the width of the 
solar peak in Fig. 1). The combination of the mea- 
sured solar and yearly modulations of the accep- 
tance may generate sidebands of equal amplitudes 
at the sidereal and anti-sidereal frequencies. 

The analysis was repeated in 3 energy ranges with 
no significant sidereal modulation observed. In the 
energy range 1 < i? < 3 EeV, with a statistics 
of almost N = 10^ events, a sidereal amplitude 
of 0.5% is measured. The expected average noise 
from the Rayleigh distribution is \/n/N ~ 0.6%, 
and an empirical estimate of this noise from the 
data at all frequencies (except for the solar band) 
gives a comparable value. 



East- West analysis 

To remove direction-independent phenomena, 
such as atmospheric and not-constant acceptance 
effects, an interesting method consists in exploit- 
ing the differences in the numbers of counts be- 
tween the East-ward and the West-ward directions 
at a given time [8]. Due to the symmetry of the 
array, the instantaneous East-ward and West-ward 
acceptances are practically equal. Assuming a har- 
monic modulation, the numbers of counts E{t) and 
W{t), from the East and West sectors respectively, 
with an average hour angle difference St with re- 
spect to the vertical, are related to the physical CR 
intensity I{t) by § ~ '^^^^rlr^- This is indepen- 
dent of the apparatus instabilities to the first order 
The first harmonic analysis of E{t) — W{t) gives 
an amplitude and a phase {rm't'D), from which 
the amplitude and phase (i?/, 0/) of the first har- 
monic of / can be derived. The background dis- 
tribution of Ri follows the Rayleigh distribution, 
given for N events by 

r N{Ri sinSt)^ - 
P{> Ri) = cxp 

Using all the data from January 2004 to March 
2007 with a declination 5 > —70°, the solar day 
modulation of 4.2% is reduced as expected, with 
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the East- West procedure, to a (0.9±0.4)% modula- 
tion. The raw sidereal amplitude is 0.6%, and after 
the East-West procedure it is (0.7 ± 0.4)%, corre- 
sponding to a chance probability of 24%. From 
this measurement we derive a 95% upper limit 
of 1.4% on the amplitude of the first harmonic 
sidereal modulation at the median primary energy 
6 X 10^'' eV. The results of the first harmonic anal- 
ysis in sidereal time with the East- West proce- 
dure, as a function of energy threshold, is shown in 
Fig. 2. AH values are compatible with fluctuations, 
the minimum Rayleigh probability being 0.6%. 

Exposure-based study of the RA distri- 
butions 

The bidimensionnal sky exposure for a given set of 
events may be obtained through various ways [9]. 
In order to account for acceptance effects without 
absorbing real anisotropies, these methods use the 
observed zenith angle distribution of the events, 
but not the azimuth distribution. A small 60°- 
periodic azimuthal modulation of the low-energy 
event rate at zenith angles > 45°, due to the ar- 
ray geometry, may however be taken into account. 
Different models for the time evolution of the ac- 
ceptance were implemented. They can use the ob- 
served event rates, the array growth and dead times 
(using the instantaneous low-level trigger activity 
of individual tanks), as well as simple parame- 
terisations of the pressure and temperature depen- 
dance of the acceptance [7]. The derived exposure 
maps reproduce well the few % RA modulation of 
the event rate which is observed when using a few 
months of data only. This modulation vanishes al- 
most completely using the full 2-year data set of 
2005-2006. A systematic uncertainty on the pre- 
dicted RA modulation of ^ 0.4% is derived from 
the differences between exposure estimates, for the 
2-year period considered. 

We describe the results of a Rayleigh analysis 
made using the 2005-2006 T5 data set with 9 < 
60°. The exposure estimate models the time evo- 
lution of the acceptance using the detector growth 
and dead time data. Fig. 3 presents a scan over 
energy of the Rayleigh amplitude, phase and sig- 
nificance. No significant modulation is detected 
throughout the scan. 



I 0.01 

0.001 L-i ' ■ — ' — 

1 10 
EthlEeV] 

350 h ' !•■-.. . ' ' ' — — 

o 300 ■ 

S§ 250 - 

^ 200 - 

f 150 - ■ _ 

& 100 - ■ . 

50 - 

' • — ' — — 

1 10 

1 ir. ■ I 



a. 



0.1 ' • ' — • • • — I 

1 10 
E,h[EeVl 

Figure 3: Rayleigh scan over energy, correcting for 
the detector's geometrical acceptance. The dotted 
fine represents the raw Rayleigh 95% probability 
corresponding to the statistics associated to each 
energy threshold. 

The first lines of table 2 summarize the observed 
modulations in two energy ranges, 1 < i;^ < 3 
EeV and 3 < < 10 EeV, as well as in two 
declination bands for the former energy range. In 
order to be conservative, upper limits are derived 
taking into account the observed residual RA mod- 
ulation in the data. Given a number N of events, 
with an exposure-corrected RA modulation Ajata^ 
one estimates from simulation, for each possible 
amplitude ^o, the distribution p^o(^) of the fitted 
amplitude with the available statistics N. The am- 
plitude Aq such that ^ dApAoiA) — 0.95 is a 
relevant 95 % upper limit to a RA modulation. 

The Rayleigh amplitude is below 1% at 1 EeV and 
an upper limit of 1.4% is derived for 1 < < 

3 EeV (there is a systematics of 0.4%i due to the 
exposure estimation, as mentionned above). 
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Energy 
[EeV] 


Declination 
range 


Number of events 
(expected noise) 


Amplitude 
(phase) 


Antisidereal 
amplitude 


95% Upper 
limit 


1<£;<3 


All 


69641 (0.7%) 


0.7% (330°) 


0.6% 


1.4% 


1< £; < 3 


[-70°, -25°] 


33470(1.0%) 


0.6% (350°) 


1.8% 


1.7% 


1< ^ < 3 


[-25°, 25°] 


31144(1.0%) 


0.9% (310°) 


0.1% 


2.0% 


3 < < 10 


All 


7722 (2.0%) 


1.2% (100°) 


1.1% 


3.2% 


E>10 


All 


1437 (4.7%) 


3.5% (70°) 


4.2% 


8.6% 



Table 1: Rayleigh analysis in different energy and declination ranges, for the periods 2005-2006 (1 < E < 
10 EeV) and 2004-March 2007 (E > 10 EeV). For each subset, we give the number N of events and the 
Rayleigh noise y/n/N ; the Rayleigh amplitude and phase; the anti-sidereal amplitude derived from Fourier 
analysis on the modified time tmod; and the 95% upper limit on a first harmonic modulation. 



Above 10 EeV, the statistics accumulated so far is 
~ 10^ events, so that %-level systematics are not 
an issue. Furthermore, as the acceptance is satu- 
rated, the use of the CIC method for energy cali- 
bration guarantees the flatness of the cos^ 6 distri- 
bution of the events [10]. Therefore one can use 
an analytic exposure estimation, given for exam- 
ple in [11], and use all the T5 events with a re- 
constructed energy above lO'^^ eV recorded since 
2004. The results of a Rayleigh analysis are dis- 
played on the last Une of table 2 : an upper limit of 
8.6% on a RA modulation can be set. 

Conclusions 

Using the first years of Auger South data, we have 
searched for large-scale patterns in the sky distri- 
bution of cosmic rays, with special emphasis on 
those above 1 EeV. A series of complementary 
analyses has been used to study the systematics in- 
duced by detector acceptance variations and insta- 
bilities. At EeV energies, the RA distribution is 
remarkably compatible with an isotropic sky, and 
an upper limit on the first harmonic modulation 
of 1.4% in the energy range 1 < E < 3 EeV is 
set. This does not confirm the 4% RA modulation 
found by the AGASA experiment (however the 
sky regions covered by both experiments are not 
the same). The results presented here concern the 
search for right ascension modulations, but stud- 
ies are also going on to search for bidimensionnal 
patterns such as a possible dipole. 
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